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TITLES 

TYPE: 

REFERENCES: 



DESCRIPTION S 
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DIGITAL COMPUTER 

REVISED LIBRARY ROOTIg E Ml^ I83 
Williams Memory Routine for Linear Programing by the 
Simplex Method (SADOI or DOl) 
Complete Floating Point Program 

"Introduction to Linear Programming", Charnes, Cooper, 
and Henderson; Wiley (1953). 

"Activity Analysis of Production and Allocation", Koopmans, 
T. C. editor; Covles Commission Monograph No. 13 , Wiley 

(1951). 

"Linear Programming Methods", Heapfcy and Candler; Iowa 

State College Press (1958). 

"Linear Programming" Kurt Eissman, Quarterly of Applied 

Mathematics , Vol. 13, October, 1955- 

This program is suitable for maximization or minimization 

problems. 

Maximization Problem 

Given a linear functional: 

(1) Z = C^ + CgXg + ... + C^ + ... + C n X n 
which is to be maximized subject to the restraints 

A u X 1 + A 12 X 2+ ... + A ln X n <B 1 

. *. » * * •• • • • ■•••• 

(2) A^l + V^ — +*^<\ 
A ml X l + A m2 X 2 + — +A «A|S B . 

and X ± > 0, 

Whenever an inequality appears a new X, is added to 
that relation making it an equality. The C ± value 
for these new %*b are equal to zero. Relation (2) 
is then written: 

Vl + •• • +- W n + X n+1 + ' ° + — + P - h 

.«••••*•« • ••••• • 

...... ..* •••«.. • 

(3) A^ + ... + A 2n x n + + X^i ... + O - B g 

\i x i + ••' V» + +-... + + ^ - B m 
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Relation (l) is unchanged, i.e. 

z - c A ♦ ... + c n x n + o x n+1 + o x n+2 + ... + o x p 

Consider the columns of coefficients of the X ' a. as 
matrix vectors % then (3) can be written: 



f ) ;X X 



"" ""* 




r - 








»• - 




* «" 


^ 


*»>- 


A ll 




^2 




Am 




i 









B l 


*21 




A 22 




*2n 









• 




• 


• 


+ X 2 


• 


T • . .X 

n 


• 


+X ^ 

n+1 


• 


+...+X 

p 


• 


= 


• 


« 




• 




• 




r - 









• 







J" 2 




A 
mn 









1 




■B 
m 
- 



Designating these vectors by P we can write relation 

(5) X,P, + X Q P + ... + x P = P„ 
X i 2 2, p p o 

wfcere P Q is the vector of the coefficient, appearing on 

the right side of the equations (3). The elements of 

f> should b# non-negative. Note that the left side of 

the equations have been rearranged so that the last » 

of the P^'s are unit vectors and form the so-called 

"Unit Basis". They are the form. 



1 


(6) 




1 







1 







Together they can be thought of as forming an identity 
matrix of order ms 



(7) 




1 




. 
. 
. 



. . . ^ 
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If a complete set of unit vectors cannot be found in 
(5)j e«g. if one of the relations in (2) is an equality 
and requires no new X. (j > n) to make it an equality, 
then the appropriate unit vectors are introduced to 
complete the unit "basis and are referred to as artificial 
vectors • 

Set the C. value for these artificial vectors equal to 
-M, where M is an extremely large number, which will 
overpower all other nmmbers involved. For example, 
one might use 9 x l(r if all the other numbers were 
of the order of magnitude of less than 1000. The 
artificial vectors are a mathematical device. only and 
should have no significance in the problem- • 

Now relation (l) becomes: 

(8) Z = C^ + ... + C n X n -MX£ - MX£ +1 - ... ->Kg +t 

where Xp , . . . , X p . are the artif iclal variables 
introduced above. Hence it is intuitive that in order 

to maximize Z, X n = X/ +1 = ••• = x £ +t = ° and they 
have no net effect on relation (3). Thift can also be 
proved in theory. 

Minimization Problem - 

Given a linear functional: 

(i) z = c^ + c 2 x 2 + ... c n x n 

which is to be minimized subject to the restraints z 
A 11 X 1 + A 12 X 2+ ... + A ln X n >B 1 



A 21 X 1 + A 22 X 2+ ..-+A 2i X n >B 2 



(9) * 



A ,X, + A ^x, + . . . + A m x > :b 
mil m2 ^ mnn— m 

Wherever an inequality appears a. new X. is subtracted 

from that relation to make it an equality. Relations 

(9) can thus be written: 



- k - 



A 11 X 1 + •- + A ln X n * X n + 1 + ° + •" + ° = B ! 



(10) 



Wl + ••• + A 2n X n +l ° " X n+2 + ° + -" +0 = B * 



2 



AX + . . . + A m X +0+...+0+X = B 
mi i mn n p m 

(p < m+n) 
C i = for all new X 1 introduced to reduce an inequality 
to an equality. 

As in the maximization problem, select a set of unit 

vectors, (the -l's introduced above do not constitute 

unit vectors), introducing sufficient artificial vectors 

to make a complete unit basis. The C.'s for the artificial 

J 
vectors in this case are +M where M corresponds to the 

above. The artificial vectors have no significance 

in this case either. 

Relation (l) now is 

(11) Z = C^ + ... + C n X n + MX^ + MX£ +1 + ... + MX/ +t 

and Xj2_ = x +1 = ••• = x / + t = ° "by a similar argument 
as before. 
A BRIEF DESCRIPTION OF THE SIMPLEX METHODS 

This is an iterative method which proceeds by introducing 
a non -basis vector into the basis and removing one of the 
old basis vectors, each time bringing the functional closer 
to its optimum value. A best solution is reached in a 
finite number of steps. 
Define the elements of P Q as A^, A g0 , A^, ..., A^. Then 

if P, thru P are the basis vectors, X, = A, rt . X„ = A_, 
J- m * l 10' 2 20* 

X 3 = A 30' '"' X m = A m0 and a11 the other x ' s are ° so 

Z = C X + C Xp + ... + C X V = C.A.. +C A on +...+CA, 
i- 1 - £ <?■ mm 1 10 ? 20 mmO 

Then Z is defined as 

Z = C,A + C A , + ... + C A , 

J 1 lj 2 2j m m j 

where P. is a non-basis vector. 

j 
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In the simplex method it turns out that when a non-^basis 
P, is introduced into the basis, the new value of the 
functional, Z', is given by: 

z- - s z - y^ - c k ) 

Since the X, is considered to be nositive P, must be 
k k 

chosen such that 

\ - c k < ° 

to increase the value of Z in a maximization or Z. - C, > 

k k 

to decrease the value of Z in a minimization. 

When all the possible Z - C. are positive the maximum 

J J 
solution has been attained or when they are all negative 

the minimum solution. Thus the end tests for maximum 

or minimum problems are to require all positive or all 

negative (Z. - C.)'s respectively. 
J J 

For purposes of using the same routine to compute the 
functional and the Z. - C. values, the functional Z is 
called Z with C Q = so that Z = Z - C . 



METHOD OF USING THIS PROGRAM; 

(a) Scaling 



To insure a minimum of roundoff error, the equations 

should be multiplied by a suitable power of ten, so 

that the elements of the matrix are of the same order 

of magnitude. For C. £ such scaling will also change 

j 
the values of the C.'s. . E.g. scaling by rows alters 



nothingj scaling a column K by (is equivalent to the 
following change of variable. 

X k = (^ so now Cy =^C k , A lk ' -fi^, ... A— = <f A 



mk 



Therefore; if the kth column is scaled by C then C of Z 

_ k 

must also be scaled by \ . 

(b) Preparing the problem 

Those vectors not appearing in the unit basis will be 

referred to as structural vectors. Number the structural 
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vectors in order starting with 1, 2, ... n where the 
P vector is as described in equation (5)5 number the 
basis vectors in order starting with n + 1, n + 2, . .., 
n + m. 

(c) The coefficient tape 

The routine performs its operations in a floating decimal 
representation so the P v 's and C, ' s have to be punched in 
this form. 

Each number is represented in the form A x ICr where 
1 > |A| > l/lO, and 6k > P > - 64. In a single register 
of the memory the number A is placed in the 33 most 
significant binary digits (& Q , a.^ ... sl ) in the same 
way as an ordinary fraction is placed in the entire 
register. An accuracy of between 8 and 9 decimal digits 
is therefore achieved. The exponent p is stored as 
the integer p + 6k in the J least significant digits 
of the same register. 

The only exception to the above rules is the number zero 
which cannot, of course, be represented as A x ICr with 
|A| > l/lO. For this reason zero is handled in a special 
way. It is represented as a number A = and p = -64. 
This representation happens to correspond exactly with 
the ordinary machine representation of zero. 

All data are prepared by translating into signed floating 
decimal designation. The first signed set refers to 
the coefficient, A, and the second to the exponent, p# 

E.g. 439.0 is punched as +439 +03 
.4390 . is punched as +4-39 +00 
.00439 is punched as +439 -02 
-43.90 is punched as -439 +02 

The coefficient tape is made by punching in order the 

elements of the structural vectors starting with P , 

P., , ... P . The basis vectors are not punched. 
1' n 
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Following the structural vectors , on the same tape, the 
C. values are punched starting with C Q (=0), C, , C p , .. ., 
C . C . , n , .... G , which are all the C . values corresponding 
to all the vectors, both structural and basis. 

There are no end marks required on the coefficient tape 
because the parameters, read in previous to this, specify 
the number of characters the program will read. 

(d) The parameter tape 

The parameter tape has the following form: 

005+ 

00F 00 jF 
OOF OOiF 
0019+ 
00F OOtF 
OOF OObF 
2k 999N 
where 

j = number of strucutral vectors not including 

P n or basis vectors, 
i = number of basis vector (this is also the number 

of restraints) 
t = frequency of the intermediate printout^ t can be 

0, 1, 2, ... (See f below), 
b = for a minimization problem 
= 1 for a maximization problem. 

The capacity of the routine is defined by : 

^(i +' j + 1) + ij < 712 

(e) Reading in the program 

The master tape is read in until the first black switch 
stop and then the parameter tape is read in. The rest 
of the master tape is then read in and will stop, after 
which the coefficient tape is read in. 

(f) The punch-out 

For intermediate results the X's and (Z. - C.)'s will 

J J 



RUNNING TIME: 



EXAMPLE: 
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be punched out every t iterationfas specified on the 

parameter tape. The number of the vector in sexadecimal 

notation will be punched followed by the X'S for the 

current basis and by the (Z. - C.)'s for the current 

J J 

non-basis vectbrs, in floating decimal form. 

If the problem converges, the final punch-out will consist 

first of the numbers of the final basis vectors in 

sexadecimal notation followed by the respective X values 

in floating decimal form. Then the final set of non-basis 

vectors and their respective Z. - C. values. 

J J 

If the problem diverges, the machine will stop on an FF 

order. 

Running time is difficult to predict from i and j above. 

Following are some samples (all with only final printout). 

In general, requiring the intermediate print-out of 

every iteration doubles machine time . 



Maximization 
i j time 
19 11 8 min 
18.25 20 min 
*18 25 15 min 



Minimization 
i J time 
8. 18 5 min 
15 31 1^ min 



Problems utilizing nearly the full capacity have almost 

always fallen in the range of 10 to 31 minutes (with 

only final printout). 

Consider the following problem: 

Maximize s 

Z = -5X X + 38Xg + ^X + 8X^ 

subject to the restraints: 
X. > 



* Same coefficients as line above but different C . . 
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+ 2X^ = 10 



14 



(a) X x - 5X 2 - 7X 5 + 2X k = 

X 2 + 6X 3 " ^ X 4 - k 
the last inequality would be written as: 

X£ + 6X 5 - 4x^ + X 5 = 4 

but Z will remain unchanged because CL = 0. 

5 
One would then rearrange the subscripts so that the 

unit vectors would appear on the right. Thus X p 

becomes X-^ or X 2 -*X 1 ', X -»X 2 ', X^ -> X ' , X ± -> X ' 

X 5" X 6' 
producing 

Z = 3ax i * + 4X 2 ' + 8X * - 5X ' + 0X 6 ' 



and 



+ 2X. 



1 ' -"3 
(b) -5X X ' - 7X 2 » + 2X ' + X 5 ' 

X,' + 6X * - 4x ' 
I 2 t> 

which can be written* 



(c) 



.= 10 
= 14 

+ XS = 4 



1 





2 





5- 


-7 


2 


1 


1 


6 


-4 






o" 




V 




10 







V 


= 


14 


1 




¥ 




>J 



L D J 

However, this does not give a complete unit basis 
on the rights the desired form is: 



(d) 



1 





2 


1 





-5 


-7 


2 





1 


1 


6 


-4 









[Y 



V 
V 



10 

Ik 



\1 



I 

61 
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which can he achieved by making (y = -100 (or an even 
larger negative number). 



Then 



z' = 38X3/ + ^x 2 » + 8x 5 ' - lOOX^' - 5x 5 ' 



It would appear intuitively obvious that in order to make 
Z* maximum, X^ = necessarily, because X^ > (this 
argument does not apply to C ' = -5 because it is of the 
same order of magnitude as C^, Cg 1 and C^ 1 ). 

Then if X^' =0, relations (c) and (d) are "equivalent" 
in the sense that when the matrices are multiplied out, 
X> * appears only in the top row, thus: 

X 1 ' +.2X 5 + X^' = 10 

but by making <y = -100, X^' = is forced which is 
equivalent to: 

X, ' + 2X ' = 10 
1 5 

In order to punch the data and parameter tapes we define: 



P = 



10 

Ik 

and the structural vectors are: 



p l = 


" 1 


P 2 = 





• V 


2 




-5 


t 


-7 


9 


2 




_ 1_ 




_6_ 




.-*. 



and F,,.T f and P^ are the unit vectors 



r l" 










i 


1 


, 








and 



respectively* 
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The data 


tape is 


as follows: 








+1 +0P 


| 


+38 +02 ) 






+lk +02 


r° 


+ k +01 








+ k +01 


1 


+ 8 +01 
- 1 +03 


» 


V 


then C/- 1 


+ 1 +oT 




- 5 +01 








- 5 +01 


' p i 


+ +00 








+ 1 +01 
J 




r 






+ +00 


) 








- 7 +01 


U 











+ 6 +01 



+ 2 +01 



+ 2 +01 
- k +01 



f P, 



y 



+ +00 



Note that P, thru P^ are not punched ,, as these are generated 
by the program. 

The parameter tape is : 



005+ 

00F 003F 
OOF 003F 
0019+ 
OOF 001F 
OOF 001F 
2^999N 



(3 structural vectors ) 
(3 basis vectors) 

(print after each iteration) 
(Maximize Z) 
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The output has the form 

004 +100000000 +02 

005 +1^0000000 +02 

006 +400000000 +01 

000 -107000000 +o4 

001 -113000000 +03 

002 +310000000 +02 

003 -218000000 +03 

003 +500000000 +01 

005 +400000001 +01 

006 +240000000 +02 

000 +200000000 +02 

001 -400000003 +01 

002 +310000000 +02 

004 +109000000 +03 



X, 



z i- c : 

z 2 " c ; 
z^ - c. 



X. 



z 1 -c ] 

Z 2 ~h 

Zi, - c, 



1st iteration 



2nd iteration 



003 +100000000 +01 

005 +520000000 +02 

001 +800000000 +01 

000 +520000002 +02 

006 +133333335 +01 

002 +390000000 +02 

004 +111666667 +03 



Z 6 ' C t 

z 2 - c ; 
z,. - c, 



3rd iteration and 
final solution 



U 
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As can be seen the problem required 3 iterations 

The final solution contains X, = 8, X, = 1, and 

X,_ = 52 with the maximum Z = 52 . 
5 
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PARAMETERS 

""* 

1 - — Temporary storage for X 1 and A 1 

2 _ 

3 Location of floating accumulator 
k Location of Library Routine A 1 

5 Number of j vectors, j = 1, . . . , i 

6 Number of i vectors, i = a + 1, . . . , p 
.7 Number of coefficients to be input 

8 Location of temporary storage for P. 

9 Location of C. 

10 Location of C. 

11 Location of Z x 

12 00 F 
00 1S5 

13 Location of X. . (j = 0, . . . , m) 

14 j vector designations 

15 i vector designations 

16 00 311F 
00 311P 

17 00 S6 
00 S6 

18 00 F 
00 311F 

19 Print parameter 

20 Maximization of minimization parameter 

INTERLUDE TO SET PARAMETERS 

00 21K 

50 12F 
75 6F 

1 S5 F 
Lk 12F 

2 L^ 6F 
^0 7F 

3 LO 12F 
LO 6F 

k Lk 18F 

40 9F 



5 


L4 12F 




40 10F 


6 


l4 6f 




1*0 l4F 


7 


L4 12F 




40 15F 


8 


l4 6f 




40 11F 


9 


L4 12F 




40 13F 


10 


L4 12F 




40 8F 


11 


50 12L 




26 999F 


12 


00 F 




00 L 


13 


00 F 




26 L 




26 IN 



INTERLUDE TO GENERATE VECTOR INTEGERS 






Ll 6f 




LO 5F 


1 


40 8L 




4l SF 


2 


F5 SF 




40 1SF 


3 


F5 8L 




40 8L 


4 


36 9L 




L5 2L 


5 


L4 7L 




40 2L 


6 


26 2L 




00 F 


7 


00 IF 




00 IF 



6 


00 F 




00 F 


9 


50 10L 




26 999F 


10 


00 F 




00 L 


11 


00 F 




26 L 




26 IN 



llbtes See "An Introduction to Linear Prograianiing?, Cnarnes, Cooper and Henderson, 
K' - coliam integer 



p r 



• row integer 

-■ column vector to be moved into basis 



i LOCATIOH 1 


ORDER 


NOTES PAGE 1 ] 





kl 11F 
L5 16F 






1 


k2 kL 








50 1L 




Enter Routine A 1 


2 


26 Sk 
OK SN 






3 


IX S6 

8k f 






k 


8S 3F 

15 F 






5 


.17 SK 




* Z J =Zx U C i 




8k 3F 




\ 


6 


8S 3F 
12 4L 






7 


85 3F 








00 S9 




- z j " c j 


8 . 


OS SS 








8J 69L 






9 


03 3L 
2K S5 






10 


8K F 
8S kF 






11 


21 1SS 




-j Test Z - C. (+) reject, (-) store new 


12 


80 ^F 
83 14L 








21 1SS 






13 


8S kF 








8J 71L 




Store K f 


Ik 


8J 73L 
23 HL 




Step K f 


■ 


■ 


:_!. 
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LOCATION 



ORDER 



NOTES 
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M 15 



15 


85 ^F 




82 68L 


16 - 


8J 73L 




3X S6 


17 


8K F 




3N F 


18 


83 25L 




31 F 


19 


83 25L 




83 25L 


20 


31 311F 




36 F 


21 


8S 3F 




Qk 120L 


22 


82 23L 




81 3F 


23 


83 25L 




81 3F 


2k 


8S 120L 




8J 81L 


25 


8J 73L 




33 17L 


26 


85 120L 




80 121L 


27 


82 65L 




82 65L 


28 


1»-K S6 




45 F 


29 


kS S8 




42 28L 


50 


8J 85L 




5K S6 


31 


8K F 




5S F 


32 


53 31L 




8K IF 



Is prograa over? 
(+) yes, (-) no 



Is x.„ = 0? 

1A 



(+) yes, (-) no 



Is x iK < 0? (+) yes, (-) no 



Waste 
-©. 



< 



Is © - © > 0? (+) store new 0, (-) 
— reject 



Store r 1 
Step r' 



Has problem diverged? (+) yes, (-) no 



Waste 



Store P , 

IV 



Clear column K f to zeros 



^L OCATION 
33 

35 
36 
37 
38 

39 
40 

41 
k2 

^3 
1* 

45 

*7 

49 
50 



ORDER 
82 33L 
8S F 
6K SN 

85 F 

86 20F 
6S SJ 
8J 93L 
62 5^L 
3K SN 
4K S6 
41 S8 
37 SJ 
kh F 
4S F 
43 38L 
8J 95L 
32 37L 
5K SN 
55 SJ 
8S F 
8J 97L 
53 42L 
8J 99L 
2K S6 

8f if 
8j 56L 
23 3HF 
89 9F 
22 45L 
22 45L 
8F IF 
7K SH 
8J 60L 
75 SS 
89 9F 
73 i^9L 



NOTES 
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M 15 



Waste Store 1 at (r'j K') 



• x' 



K3 



~*'lJ- X iJ- ( ^ )X : 



IK 



VK 



— Restore x' 



£J 



-Print X's 



Waste 



Print Z, - C. 



LOCATION 


ORDER 


NOTES PAGE k M 15 




51 


73 J+9L 
8j 6kL 




Waste 




52 


OF F 
50 52L 








53 


26 Sk 
4K S7 








> 


88 F 




- Input 




: 


ks 3HF 








55 


43 5^L 
8j L 








56 


92 131F 
L5 SL 








57 


00 28F 
82 12F 




-Print \ designation 




58 


92 963F 
F5 56L 








59 


ko 56L 

26 29S^ 








60 


92 131F 
L5 SF 








61 


00 28F 








62 


82 12F 
92 96 3F 
F5 60L 




"Print Z. - C. designations 
J J 




63 


ko 60L 
26 29S4 




. 




6k 


LI 19F 
kO 12F 




■Reset print counter 




65 


26 83L 

8j 66l 








66 


FF F 
00 52L 




» 




61 


L5 66L 
k2 51L 




■Prepare to stop, program 




68 


26 116L 
8j 67L 






JL 



LOCATION 



ORDER 
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69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 
84 

85 
86 



L5 4l 
L4 6f 
4o 4l 
26 29S4 
F5 11F 
40 13F 
26 29S4 
00 F 
F5 11F 

40 11F 
26 29SV 
00 F 

50 13F 
75 6F 
L5 16F 
S4 F 
42 17L 
42 18L 
42 20L 
42 28L 
42 31L 

41 11F 
26 114L 
00 F 

L5 11F 

40 18F 
26 29S4 
00 F 

41 S3 
50 83L 
26 S4 
82 16L 
L5 28L 
L4 18L 

42 33L 
42 89L 



— From 8L 



- Store K 1 



■^Step K' of r' 



— From l6L 



- Store r 1 



LOCATION 
87 

88 

89 
90 
91 
92 
93 
94 
95 
96 
97 
98 

99 

100 

101 

102 

103 

104 



ORDER 
L5 16F 
L4 18F 
42 3^L 
42 42L 
22 89L 
I£ F 
40 20F 
L5 16F 
42 39L 
46 39L 
26 29S4 
00 F 
L5 34L 
L4 6f 
40 3^L 
26 29S4 
L5 39L 
L4 17F 
40 39L 
26 29S4 
L5 42L 
L4 6f 
4o 42L 
26 29S4 
L5 9F 
L4 13F 
42 103L 
42 105L 
L5 10F 
L4 18F 
42 104L 
42 106L 
22 103L 
50 F 
22 104L 
L5 F 



NOTES 
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M 15 



.From 30L 



Store X 



rK 



» From 36L 



1 



^ From 40L 



From 43L 



Interchange i and j designation; also 
c and c 



LOCATION 
105 

106 

107 
108 
109 
110 

111 



113 
114 

115 
116 

117 
118 

119 
120 
121 



ORDER 
22 10% 
40 F 
S5 F 
40 F 
LI 55L 
40 53L 
36 L 
L5 l4F 
iA 13F 
42 103L 
42 105L 
L5 15F 
iA 18F 
42 104L 
42 106L 
L5 121L 

40 120L 
22 103L 
26 29S4 
00 F 

36 ll6L 
26 29S4 
L5 l4F 
42 60L 

L5 15F 

42 56L 

41 S3 
50 118L 
26 S4 
82 44L-. 
40 F 

00 127F 
40 F 
00 127F 
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Binary switch 



For non-zero print is modified 



The following interlude a et$ the prist and nlao aaxiBisation or minimizati 



on. 



LOCATION 



3 

k 

5 
6 

7 
8 

9 

10 
11 
12 



ORDER 



H0TB3 
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00 1^3K 
LI 19F 
32 2L 
40 12F 
L5 8L 
40 135F 
LI 20F 
36 6l 
L5 9L 
40 32F 
L5 1«L 
40 33F 
26 6L 
50 7L 
26 999F 
00 F 
00 L 
F5 12F 
40 12F 
25 IBS 
,80 4f 
83 35F 

25 1SS 
00 F 

26 L 
26 IK 

Routine A 1 
20 T5W 
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